Sub-diffraction Focussing using Quadratic Measure Eigenmodes 
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For over a century diffraction theory has been thought to limit the resolution of focusing and imag- 
ing in the optical domain. The size of the smallest spot achievable is inversely proportional to the 
range of spatial wavevectors available. Here, we show that it is possible to locally beat the diffrac- 
tion limit at the expense of efficiency. The method is based on the linearity of Maxwell's equations 
and that the interaction between light and its surroundings may be considered quadratic in nature 
with respect to the electromagnetic fields. We represent the intensity and spot size as a quadratic 
measure with associated eigenmodes. Using a dynamic diffract ive optical element, we demonstrate 
optical focussing to an area 4 times smaller than the diffraction limit. The generic method may 
be applied to numerous physical phenomena relating to linear and measurable properties of the 
electromagnetic field that can be expressed in a quadratic form. 

PACS numbers: 42. 25. Fx, 42.40.Eq 



Introduction The diffraction limit may be seen to 
originate from the Heisenberg uncertainty principle [lj|2]. 
The concept of super-resolution to beat the so termed 
"diffraction limit" is typically associated with the reten- 
tion of evanescent fields which are needed to retain all k- 
space components to realise "perfect imaging" [3ji^ . The 
concept of a superlens [6 is typically based on the recov- 
ery of evanescent information by judicious use of suitable 
metamaterials - artificial materials exhibiting a negative 
refractive index [7]. However, it would be advantageous 
to explore methods of overcoming the diffraction limit 
in the far field, even if over a restricted range, as this 
is amenable to the vast majority of optical systems. In- 
triguingly, recent work on band-limited functions oscillat- 
ing faster than the highest Fourier components of which 
they are composed, so termed superoscillations, shows 
that there may be routes to sub-diffraction imaging in the 
optical far-field without the need to retain rapidly decay- 
ing evanescent waves [8HTT1. However, with superoscilla- 
tions, we achieve sub-diffraction features at the expense 
of having most of the energy contained in low frequency 
Fourier features (side-bands) that have many orders of 
magnitude higher amplitude than any sub- diffract ive fea- 
ture we may wish to utilise. Speckled fields, originating 
from strong (diff user- like) phase aberrations, can also de- 
liver random positioned sub-diffraction patterns or hot 
spots. Their usefulness is limited due to their random- 
ness and the lack of control of the distribution of energy 
between the speckles. It thus remains an open question 
as to how to potentially realise sub-diffraction focusing 
and imaging in a generic manner that is applicable as 
this would open up a wealth of new scientific directions. 

Here we address the problem of realising the intensity 
of an optical field such that we produce a focussed spot 
smaller that the diffraction limit dictated by the optical 
system employed. The most intuitive way to describe our 
approach is as a superposition of incident fields that is op- 
timised to achieve this smallest spot. Naturally, multiple 



techniques can be employed for this optimisation process 
ranging from a genetic algorithm to a steepest descent 
method. The major challenge encountered in any such 
optimisation and engineering of optical intensities is the 
fact that electromagnetic waves interfere. Indeed, the 
intensity profile arising from the superposition of mul- 
tiple fields is a complex interference pattern depending 
not only upon the intensities of the considered waves but 
also upon their relative phase fronts. This interference 
pattern not only makes the search for an optimum beam 
problematic but crucially renders the superposition found 
unreliable as the different algorithms may converge on 
local minima which are unstable with respect to the dif- 
ferent initial parameters in the problem. In contrast, our 
proposed method yields a unique solution to the prob- 
lem and directly yields sub-diffraction optical features in 
the far field of our designated optical system. As might 
be expected, our method gives insights into the area of 
superoscillations which we shall return to in the conclu- 
sions. 

Our method is based on two fundamental properties 
of the electromagnetic field and its interactions. Firstly, 
the approach relies on the linearity of the electromag- 
netic fields i.e. that the sum of two solutions of Maxwell's 
equations is itself a solution of them. As we consider free 
space propagation, this criteria is satisfied. The second 
property relates to the interaction of the electromagnetic 
field with its environment. All such interactions can be 
written in the form of quadratic expressions with respect 
to the electric and magnetic fields. Examples include the 
energy density, the energy fiow and Maxwell's stress ten- 
sor. This allows us to designate appropriate "quadratic 
measure" eigenmodes to various parameters {e.g. spot 
size) and subsequently ascertain the minimum eigenvalue 
which, in the case of a spot size operator yields a sub- 
diffraction optical feature. We first present the theory 
underlying our approach and then demonstrate experi- 
mentally the applicability of our method to achieve sub- 
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diffraction focussing. 

Quadratic measure eigenmodes (QME) Our method 
assumes monochromatic solutions of free space Maxweh's 
equations, with E and H the electric and magnetic vector 
fields and with eo and /io the vacuum permittivity and 
permeability. These solutions can be written in an inte- 
gral form linking the electromagnetic fields on the surface 
A with the fields at any position r 

Fu(v)= [ P,,(r,rO^.(r0^5' (1) 

J A 

where a/2J^ = (y^E, yT^H) is a shorthand for the two 
electromagnetic fields having six scalar components. 
The integration kernel P^^ corresponds to a propaga- 
tion operator giving rise to different vector diffraction 
integrals such as Huygens, Kirchhoff and Stratton-Chu. 
Crucially all linear and measurable properties of the elec- 
tromagnetic field can be expressed as quadratic forms of 
the local vector fields and are therefore termed quadratic 
measures. For instance, the time averaged energy density 
of the field is proportional to J^* • J* = 1/2 (eoE* • E + 
/ioH* -H) while the energy flux to 1/2(E* x H + E x H*). 
The asterisk * stands for the complex conjugate. Inte- 
grating the first quantity over a volume determines the 
total electromagnetic energy in this volume while inte- 
grating the normal energy flux across a surface the in- 
tensity of the light field incident on this surface. All the 
quadratic measures can be represented in a compact way 
by considering the integral 

M,= [ T^KuvJ'vdr = {T\k\T)v (2) 
Jv 

where the kernel Kuv = i^lu Hermitian where f the 
adjoint operator including boundary effects for finite vol- 
umes. Table |T] enumerates some operators associated to 
common quadratic measures. The integrand part of all 
these quadratic measures correspond to the conserving 
densities which together with the associated currents is 
Lorentz invariant ^2 - The volume, over which the inte- 
gral is taken, does not need to be the whole space and 
can be a region of space, a surface, a curve or simply 
multiple points. To account for this general integra- 
tion volume, we broadly term it the region of interest 
(ROI) in the following. Finally, using the general defi- 
nition of the quadratic measure it is possible to define a 
Hilbert sub-space, over the solutions of Maxwell's equa- 
tions, with the energy operator (EO) defining the inner 
product. Further, any general quadratic measure defined 
by ([2| can be represented in this Hilbert space by means 
of its spectrum of eigenvalues and eigenfunctions defined 
by XTu = i^uvJ^v Depending on the operator hZuv^ the 
eigenvalues A form a continuous or discrete real valued 
spectra which can be ordered. This gives direct access 
to the solution of Maxwell's equations with the largest 
or smallest measure. The eigenfunctions are orthogo- 
nal to each other ensuring simultaneous linearity in both 



field and measure. In the following, we study the case of 
different quadratic measure operators and their spectral 
decomposition into modes which we term the quadratic 
measure eigenmodes ( QME). The convention for operator 
labelling we adopt is to use the shorthand QME followed 
by a colon and a shorthand of the operator name. 



Operator 


2»7^ iiU'uv'^ V 


EO 


eoE* . E + /xoH* • H 


10 


(E* X H + E X H*) • Ufc 


sso 


r2(E* X H + E X H*) -Ufe 


LMO 


eoE*.(iafe)E + /ioH*-(iafc)H 


OAMO 


eoE* • {ir x V)fcE + /ioH* • {ir x V)fcH 


CSO 


i(E* H - H* E) 



TABLE I: Common quadratic measure operators including 
the energy operator (EO), intensity operator (10), spot size 
operator (SSO), linear momentum operator (LMO), orbital 
angular momentum operator (OAMO) and circular spin op- 
erator (CSO). The vector operators include the subscript k 
indicating the different coordinates and the associated unit 
vectors. 

Practical example: smallest focal spot In our practi- 
cal example we utilised the intensity operator QME:IO 
and the spot size operator QME: SSO defined in the fol- 
lowing to engineer the size of a laser focus. The QME:IO 
measures the electromagnetic energy flow across a surface 
A: 

^(0) ^ \ j (E* X H + E X H*) • nd5 (3) 

where n is normal to the surface of interest. The eigen- 
vector decomposition of this operator can be used for 
example to maximise the optical throughput through a 
pinhole or to minimise the intensity in dark spots. Con- 
sidering a closed surface surrounding an absorbing par- 
ticle, the QME of the 10 give access to the field that 
either maximises or minimises the absorption of this par- 
ticle. The definition of the QMEiSSO is based on the 
concept of determining the spot size of a laser beam by 
measuring, keeping the total intensity constant, the sec- 
ond order momentum of its intensity distribution. 

m(2) = ]- j |r - rop(E* x H + E x H*) • nd5 (4) 

where r is the position vector and tq the centre of the 
beam. The eigenvalues of this operator measure the 
spread of the beam with respect to its centre and the 
smallest eigenvalue defines the smallest spot achievable 
in the ROI. In the following, we define a = VmJ^ as the 
width of this spot. 

For the experimental determination of both the 
QME:IO and QME: SSO, we considered a pair of an initial 
and a target plane located at the propagation distances 
z = zj and z = zt and connected through a linear optical 
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system. Crucially, a superposition of fields E{x,y^zi) = 
'^u=i ^uEu{x^ zi) {au G C, Nu G N) in the initial plane 
is rendered into a superposition of the respective propa- 
gated fields E{x,y,ZT) = Yl^li^uEu{x,y, z^y) charac- 
terised by the same set of coefficients au due to lin- 
earity of the optical system. Based on this superposi- 
tion approach and the QME:IO as defined in Eq. ([3|, 
the intensity in the target plane can be represented as 
m^°) = a*M^^)a. M^^) is N x N matrix with the ele- 
ments given by the overlap integrals 

^i'^ = / K{^. zt)Eu{x, y, ZT)dS. (5) 

^ROI 

This matrix is equivalent to the QMEiIO on the Hilbert 
subspace defined by the fields Eu{x^y^ zt)- M*^^^ is Her- 
mitian and positive-definite which implies that its eigen- 
values X^^^ {k — 1 . . . Nu) are real and positive and the 
eigenvectors v[.^^ are mutually orthogonal. Accordingly 
the largest eigenvalue Amax = max(A[.^^) and the as- 
sociated eigenvector Vmax will deliver the superposition 
E^^^[x,y,z) = Y.^l^v^^l^^uEu{x,y,z) {z = zi and z = 
zt due to linearity) which maximizes the intensity within 
the ROI. Similar to the QMEiIO, the QME:SSO as de- 
fined in Eq. Q can be written as m^^^ = b*M^^^b where 
M^^^ must be represented in the intensity normalised 
base E,ix,y,ZT) = Y.u=Mkl/>^T)Eu{x,y,z^). M(2) 
IS di N X N matrix with the elements given by 

m£)= / \v-v^\^El{x,y,z^)Eu{x,y,z^)dS. (6) 

J A 

We denote the eigenvalues of M^^^ as A^^^ and the 

eigenvectors as \ The eigenvector associated with 
the smallest eigenvalue corresponds to the smallest spot 
achievable within the ROI through the linear superposi- 
tion of the Nu fields Eu{x^y^zi) considered initially. 

Experiment Our experiments are based on the ex- 
pressions ([5| and (|6| which allowed us to determine the 
superposition coefficients for the smallest spot from a 
set of test electric fields Eu{x^y, zt) measured in the 
target plane. We used an expanded HeNe laser beam 
(P = 4 mW, A = 633 nm) to illuminate the chip of 
a phase-only spatial light modulator (SLM, type Hama- 
matsu LCOS X10468-06, 800 pixel x 600 pixel) operat- 
ing in the standard first order configuration Note 
that additional amplitude modulation can be achieved 
in a straightforward manner but is beyond the scope 
of our proof-of-principle study. Without lack of gener- 
ality we have chosen the Zernike polynomials Z^{x^y) 
{m^n G N) [14 to modulate the beam phase that is 
our test fields exhibited a phase behavior according to 
Eu{x,y^zi) (X exp(z-Z^(x, y)) where the index u enumer- 
ates the different combinations {n^m). The modulated 
beam was subsequently propagated through a spherical 
lens (focal width / = 1 m), the linear optical system, and 



then detected with a CCD camera (Easier pilot piA640- 
210gm). Since the CCD camera only detected intensi- 
ties we applied the well-known lock-in technique to the 
optical domain as described in detail in the Supplemen- 
tary information. In brief, a reference Gaussian beam, 
whose phase was oscillated in time using the SLM, was 
interfered with the test field in the target plane in or- 
der to determine both amplitude Au{x^y, zt) and phase 
(j)u{x^y^ Zt) of the test field in the target plane. The 
respective reference field parameters Aji{x^y^ zt) and 
(j)R{x^y^ z^) were independently determined using self- 
interference and an approximate evaluation of the phase 
gradient field including subsequent numerical integra- 
tion. Both the QMEiIO and the QMEiSSO were finally 
constructed from the measured parameters Au{x^y^ZT) 
and (j)u{x^y^ Zt) according to Eqs. ([5| and ([6|. Dedi- 
cated Labview and Matlab software allowed us to record 
a set of Nu test fields Eu{x^y, zt) (typically Nu = 
231 corresponding to the Zernike polynomials up to 
order n = 20) at a rate of 50 Hertz. Each test 
field required a 48 point temporal phase scan. Nu- 
merical evaluation of the QMEiIO and QME:SSO fi- 
nally delivered the required superposition E{x,y^zi) = 

A(x,^,zi)e^^(-'^'^^) = Et^i^miU^n(^,^,^i) which was 
encoded onto the SLM. The final superposition required 
simultaneous modulation of both amplitude and phase 
of the laser beam incident onto the SLM which we have 
encoded to our phase-only SLM using the approximation 
^tA(f) (described in detail elsewhere [15]. Cru- 
cially, the QMEiSSO was determined for decreasing size 
of the target ROI which allowed us to squeeze the laser 
spot size below the diffraction limit as shown in the fol- 
lowing. 

Results and discussion Figure [l] shows a set of in- 
tensity profiles I{x,y^ zt) as obtained after encoding the 
final superposition of test fields for different target ROI 
sizes. The ROI side length aRoi is indicated in the pro- 
file's left top corner in units of the Airy disk size (JAiry-disk 
which was (lAiry-disk = 63 |LLm given the laser wavelength 
A and the numerical aperture NA = 0.005 of our optical 
setup. Crucially, the intensity profiles reveal a central 
spot whose size is decreasing when the ROI size is re- 
duced. This is balanced by a redistribution of intensity 
into the area outside of the ROI. Interestingly, the ap- 
plied procedure not only aims to achieve the smallest 
spot size possible for a given set of test fields but also 
clearly aims to keep the redistributed intensity entirely 
outside the square shaped ROI. The redistributed inten- 
sity starts to evolve at cr/dAiry-disk ^ 6 (data not shown) 
and becomes predominant for cr/dAiry-disk < 4. 

We have also performed a quantitative analysis of the 
qualitative profiles shown in Fig. [l] which was based on 
the following parameters: 1. The central spot size a (de- 
termined as the FWHM of a Gaussian the central spot 
is fitted to) in units of (lAiry-disk, 2. the peak intensity 
^max of the central spot relative to the peak intensity 
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FIG. 1: 2D intensity profiles recorded in the target plane 
for decreasing ROI size. The lateral size of the profile plots is 
1.04 mm X 1.04 mm. The ROI is indicated by black and white 
dashed square. Numbers in the left top corner of the profiles 
indicate the ROI size in units of the Airy disk size ctroi = 
63 ^im. Blue to red color indicates low to high intensity. 
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FIG. 2: Quantitative analysis of the intensity profiles shown 
in Fig. [l] as a function of the width aRoi of the ROI. (a) 
Width of the central spot a as defined by Q . (b) Normalised 
maximum intensity /max- (c) Contrast ratio between the max- 
imum intensity /max,ext-ROi outside of the ROI and maximum 
intensity /max of the central spot. All lengths are measured 
in units of fJAiry-disk = 63 ^im. 



^max, Airy-disk of the Airy disk and 3. the peak inten- 
sity /ext-ROi outside of the ROI relative to the central 
spot peak intensity /max- The respective parameters were 
plotted versus aRoi/cTAiry-disk as shown in Fig. [2j^a)-(c). 
The plotting range extends to aRoi /cTAiry-disk ^ 2 where 
the measured data became very noisy since the limits 
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FIG. 3: 2D intensity profiles obtained by numerically super- 
imposing the measured fields Eu{x,y, zt) in the target plane 
as opposed to the profiles shown in Fig. ^ obtained by SLM 
encoding the final superposition in the initial plane. Figure 
specifications are equivalent to Fig. [l] 



in terms of both SLM encoding and detector sensitiv- 
ity were reached. All three parameters exhibit distinct 
changes when cr/dAiry-disk ^5: 1. The relative spot 
size cr/cr Airy- disk fastly drops to a value below 0.5 after 
having persistently decreased to approximately 0.7 (see 
Fig. [2|a)). 2. The central spot relative peak intensity 
^max/^max, Airy- disk manifests a sharp decrease of the rel- 
ative eigenspot peak intensity (see Fig. [2|b)). 3. The 
redistributed relative intensity /max,ext-ROi/^max simul- 
taneously experiences a vast increase (see Fig. [2|c)). 

Overall, our results manifest remarkable similarities to 
the concept of superoscillations in band limited functions 
[8 . The central spot size is decreased below the diffrac- 
tion limit at the expense of the spot intensity which is 
redistributed to the so-called side bands around the spot. 
It is important to remark upon the difference and com- 
monality between superoscillation and quadratic measure 
eigenmodes. Superoscillations relates to generation of 
spectrally bandwidth limited fields that produces oscil- 
lations with frequencies outside the limiting bandwidth. 
As such, super oscillating fields are based on the Fourier 
relationship between reciprocal space and real space. In 
contrast, our quadratic measure eigenmodes approach is 
a generic method that does not rely on Fourier relation- 
ship but can be used to optimise any general quadratic 
measure operator. The application of QME corresponds 
to minimising the spot size of a superposition of band- 
width limited fields and in this specific case, the QME is 
equivalent to super oscillation. 

Although we are currently limited to minimum spot 
sizes of cTAiry-disk ^ 0.5 duc to phasc-ouly SLM encod- 
ing, limited SLM resolution, and limited detector res- 
olution and sensitivity, our concept offers a huge po- 
tential. To demonstrate this we have performed a nu- 
merical superposition of the measured fields Eu{x^ zt) 
in the image plane. The resulting intensity distribution 
Iu{x,y,ZT) cx \ Y.u^u{x,y,ZT)\^ is shown in Fig. [s] and 
would be achieved with an advanced experimental config- 
uration in particular featuring high resolution amplitude 
and phase SLM encoding to realize nearly-perfect linear- 
ity in our optical system. The central spot in the right 
graph has a size of (TAiry-disk ~ 0.1 that is the Airy disk 
is beaten by an order of magnitude. 

In this letter, we have experimentally used the QME 
approach to locally generate sub-diffraction light spots 
and we demonstrate beating the diffraction limit by a 
factor of 4 in area. The theory that we employ is rig- 
orous and based on considering the spot size operator 
as a quadratic measure originating from Maxwell's equa- 
tions. Excitingly, we can define other quadratic measure 
operators to which our approach is applicable (see Ta- 
ble |l| . The generic nature of our approach means that 
it may be applied to optimise the size and contrast of 
optical dark vortices, the Raman scattering or fluores- 
cence of any samples, the optical dipole force and the 
angular/linear momentum transfer in optical manipula- 
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tion. Our approach is applicable to all linear physical 
phenomena where generalised fields interfere to give rise 
to quadratic measures. 
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